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Abstract. We introduce a discrete model for binary spin-orbit-coupled (SOC) Bose- 
Einstein condensates (BEC) trapped in a deep one-dimensional optical lattice. Two 
different types of the couplings are considered, with spatial derivatives acting inside 
each species, or between the species. The discrete system with inter-site couplings 
dominated by the SOC, while the usual hopping is negligible, emulates condensates 
composed of extremely heavy atoms, as well as those with opposite signs of the effective 
atomic masses in the two components. Stable localized composite states of miscible 
and immiscible types are constructed. The effect of the SOC on the immiscibility- 
miscibility transition in the localized complexes, which emulates the phase transition 
between insulating and conducting states in semiconductors, is studied. 
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1. Introduction 

Recently, the use of ultracold quantum gases - atomic Bose-Einstein condensates (BEC) 
and fermion gases alike - for simulating fundamental effects originating in condensed- 
matter physics has drawn much interest [1]. One of these effects is the spin-orbit 
coupling (SOC) which links the electron spin to its motion in semiconductors. The 
SOC of the Dresselhaus [2] and Rashba [3] types plays a major role in many phenomena 
and applications, including spin and anomalous Hall effects [4], topological insulators 
[5], spintronics [6, 7], spin-based quantum computations [8], etc. In contrast to the 
complex situations found in solids, the “synthetic” SOC, induced by appropriate laser 
illumination of atomic gases in the combination with a magnetic field, can be precisely 
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controlled in the experiment [9, 10, 11, 12]. Furthermore, unlike the electron spin, the 
pseudo-spin of laser-dressed atoms is not constrained by fundamental symmetries, which 
gives rise to a variety of settings imavailable in solids, such as new exotic superfluids 
[13]. The theoretical description of SOC effects in BEC is also much simpler than in 
solids [9, 14, 15]. Therefore, studies of the SOC in BEC, as well as in fermionic gases 
[16, 17], have become a vast research area, see Ref. [18] for a brief review. 

SOC in solid state is manifested by electrons moving in static electric fields 
[19, 20, 21]. The coupling results from the Zeeman interaction between the magnetic 
moment of the electron, aligned with its spin, and the magnetic field appearing in the 
reference frame moving along with the electron. To link this setting to the atomic 
gas, the electron’s spinor wave function is mapped into a pseudo-spinor mean-held 
wave function of the binary BEC in 8 'Rb, which contains atoms in two different states 
“dressed” by appropriate laser fields. In particular, making use of the electronic ground- 
state 5Si/2,F = 1 and F — 2 hyperhne manifolds, one can start with the set of four 
states, |.F, mj?) = |2, 0), 11, —1), 11, 0), 12,1), which are cyclically coupled by four strong 
fields directed along diagonals in the (x, y ) plane, with the magnetic held applied along 
z. As shown in [9], by means of an appropriate unitary transformation, which dresses the 
original atomic states, this setting may be described by a combination of the Rashba and 
Dresselhaus Hamiltonians, with independently adjustable coefficients, a and /3, acting 
on a set of two dressed states, that emulate the spin-up and spin-down polarizations of 
electrons in the semiconductor: 

H x =a ( cr x p y - (Typ x ) + /3 (er x p y + a y p x ), (1) 

where a XtV are the Pauli matrices, and p x , y are operators of the respective momentum 
components. 

In the seminal experimental work [10], a different setting was implemented in the 
condensate of 8 'Rb: a pair of Raman laser beams with wavelength A, shone along 
the same diagonal directions in the (x, y) plane and with the same vertical orientation 
of the magnetic held as mentioned above, were used to couple two states, \p> + ) = 
|F = 1, m F = 0) and l^ - ) = \F — 1, m F = —1) (the third one, | F = 1, m F = +1), which 
belongs to the same manifold, is far detuned from this pseudo-spinor set). The original 
Hamiltonian of the system includes the kinetic energy of atoms, Zeeman terms, and inter¬ 
component coupling terms modulated by factors exp (±47rix/A). The transformation of 
this Hamiltonian by unitary matrix U = exp ( 27ria z x/X ) makes it possible to eliminate 
the explicit x-dependence, simultaneously generating the effective SOC term, 

H 2 = 7°zPx, ( 2 ) 

due to the noncommutivity of U and term p 2 x / {2m) in the kinetic-energy operator. 
Further, the pseudo-spin rotation can make H 2 equivalent to the combination of the 
Rashba and Dresselhaus couplings with equal strength, i.e., a particular case of the 
generic Hamiltonian (1). 

The interplay of the SOC with intrinsic collisional nonlinearity of the pseudo-spinor 
BEC was recently considered too, in many theoretical works. In particular, it was 
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demonstrated the new ground-state phases (stripes, phase separation, etc.) can be 
created in a such nonlinear two-component systems [22], as well as tricritical points 
[23], solitons of different types [24, 25, 26, 27], including two-dimensional solitons with 
embedded vorticity [28, 29, 30], and vortex lattices [31]. New topological excitations 
emerge in spin-orbit-coupled fermionic gases too [17]. SOC gives rise to other noteworthy 
phenomena when combined with an optical-lattice (OL) potential: flattening of the 
Bloch potential [32], atomic Zitterbewegung [33], and new topological phases [34], 

Our objective here is to derive a discrete form of the SOC-BEC model for the 
pseudo-spinor condensate trapped in a deep OL potential, and then construct discrete 
solitons specific to this model. Especially interesting is the regime in which the inter-site 
coupling is dominated by SOC, while the usual hopping may be neglected. It gives rise 
to a previously unknown discrete model that makes it possible to emulate the SOC in 
BEC composed of infinitely heavy atoms, as well as emulate binary condensates with 
opposite signs of effective atomic masses (dynamics of trapped SOC BEC featuring 
an infinite effective mass in some special cases was recently addressed in non-discrete 
settings both experimentally [12] and theoretically [15]). 

We also investigate the immiscibility-miscibility (IM) transition induced by the 
SOC in localized pseudo-spin complexes. In this connection, it is relevant to mention 
that, in the absence of linear couplings between the components, the hyperhne states 
in 87 Rb are immiscible, although being close to the miscibility threshold [35, 36, 37]. 
Previously elaborated proposals to achieve a tunable immiscibility-miscibility transition 
relied on the adjustment of the scattering length via the Feshbach resonance, controlled 
by magnetic [38, 39] or optical [40, 41] fields, or on linear coupling between the different 
atomic species imposed by a radio-frequency held [42, 43]. The SOC imposes its specific 
linear mixing via spatial gradients of the wave functions, thus offering a new way to 
induce the immiscibility-miscibility transition. 

The SOC model, written in terms of coupled discrete Gross-Pitaevskii equations 
(GPEs), is introduced in Section II. Results for discrete bright-soliton complexes in this 
system and their stability, as well as for the effect of the SOC on the immiscibility- 
miscibility transition in the binary condensate, are presented in Section III. The paper 
is concluded by Section IV. 

2. The model 

2.1. Basic equations 

We consider the BEC composed of atoms in different hyperhne ground states 
corresponding, as said above, to F = 1 [10] or F — 1,2 [9]. The states are coupled 
by two [10] or four [9] Raman laser beams shone at angle 7t/4 in he x — y plane with 
respect to the system’s axis x, with the dc magnetic held applied along z. These settings 
induce the SOC in the basis of two dressed states, which emulate the spin up and down 
components. In the presence of the OL potential, Vol ; the mean-held dynamics of the 
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BEC is governed by the GPE written in the spinor form: 


lh lTt = 


2m 


+ Vol + H soc + H, 


*mix/spl 


EH, 


int 


A>. 


( 3 ) 


where T = (0 + ,0 _ ) T is the normalized spinor wave function in the dressed-state 
representation, and p 2 / (2m) represents the kinetic energy of the quasi-ID condensate. 

According to the above discussion, the Hamiltonian term (2), which accounts for 
the SOC in the setting similar to that which was experimentally realized in [10], amounts 
to 


Hsoc = 


tv 2 d 

-l — K(T z — , 

m ox 


( 4 ) 


where k is the SOC strength. In this case, we assume that the condensate is subject to 
tight confinement in the y and z directions, hence p 2 in Eq. (3) may be replaced by p 2 . 
On the other hand, starting from the more general SOC Hamiltonian (1), and assuming 
the tight confinement in the (z,x) plane, i.e., keeping only p y in p 2 in Eq. (3) and in 
(1), the SOC Hamiltonian amounts to 


TT h2 9 

H soc — ~ l — K<j x -pr, 
m oy 


( 5 ) 


with respective strength k. 

The Hamiltonian term (4) implies that the spatial derivatives accounting for the 
SOC act inside of each dressed state (we call this system an intra-SOC one), while the 
term (5) couples different dressed states by the derivatives (to be called the inter-SOC 
system). In the former and latter cases, the unconfined coordinate is denoted as x 
and y, respectively. The commutation of the SOC terms in the Hamiltonian with the 
respective coordinate gives rise to what may be considered as anomalous velocity, and 
effects produced by the SOC terms, such as the immiscibility-miscibility transition (see 
below) may be accordingly understood as mixing caused by this velocity. 

Further, the linear-mixing or splitting terms of the Hamiltonian, which are induced 
by the magnetic held, are represented, severally, by H mix = hfla x and H sp \ = Mla z 
in the intra- and inter-SOC systems, where 0 is, respectively, the Rabi coupling or 
Zeeman splitting between the two wave-function components. Lastly, H int accounts for 
the collisional intra- and inter-component interactions in the BEC. 

The fragmentation of the condensate in the deep OL potential into droplets coupled 
by tunnelling across barriers separating local potential wells leads to the replacement 
of the continuous GPEs by their discrete counterparts, which can be performed in 
essentially the same way (using the tight-binding approximation) as it was elaborated 
for the single-component BEC [44, 50]. Thus, adopting the units where h — 1, m = 1, 
we arrive at a system of two discrete GPEs for the (pseudo-) spinor wave function, 
A> n = (-0+, 0“) r , where n is the discrete coordinate, which replaces x in Eq. (4), or y in 
Eq. (5). The first system corresponds to the intra-SOC, cf. Eq. (4), with the first-order 
finite-difference derivatives acting separately on each component: 

Mn 


— n(nU+ L „/,+ a 


.(„lA 






2 





Composite localized modes in discretized spin-orbit-coupled Bose-Einstein condensates 5 


i-Qf = - C(l /} n+ 1 + W—l) + ^(ipn+l - 0„-i) + Wn\ 2 + 7 Hn ~ ^n(6) 

The other model corresponds to the inter-SOC, in which the finite-difference derivatives 
mix the components, cf. Eq. (5): 

*= - C'^n+l + V’n-l) + - C-l) + (7l|^ + | 2 + /3|d")^n ~ , 

= - G ^n +1 + VV-l) + ^(V'n+l “ V’n-l) + Wn | 2 + 7^0 ?)i>n + ^r£ 7 ‘) 

In either system n is the rescaled SOC strength, 7 , 7 ! < 0 and /3 < 0 are, respectively, 
coefficients of the intra- and inter-component interactions, which we assume here to 
be attractive, and O is the Rabi-coupling or Zeeman-splitting frequency in Eqs. ( 6 ) 
and (7), respectively. The nonlinearity coefficients are normalized by setting 7 x = —1. 

The remaining free parameters are n,/3, 7 , and fl, along with the inter-site hopping 
coefficient, C > 0. 

It is well known that the discrete equation with repulsive on-site nonlinearity can 
be transformed into its counterpart with the attractive interaction by means of the 
staggering transformation followed by the complex conjugation: 

^ = (-1 )>;, ( 8 ) 

where the asterisk stands for the complex conjugate. If the underlying nonlinearity is 
repulsive, the transformation will invert the original signs of 7 x, 7 and /?, while the sign 
of the SOC coefficient, k, will not change. 

Each system conserves two quantities, that are written here for infinite lattices: the 
norm, 

(wtr+h/’nD, ( 9 ) 

n 

and Hamiltonian, which is 

H = i~ C WViH-l + V’n^n+l) - iK On*^n +1 ~ ^n*^n+ 1 ) 

n 

+ \ (71 1^1 4 + 7 l^n | 4 ) + ^ \*l>n | 2 kn | 2 - ^*00 | + C.c. (10) 
for the intra-SOC system ( 6 ), and 

H = J2 i~ C i^n^n+l + 1 J n*i>n+ 1) + ^n* (V’n+l * V’n-l) 
n 

+ \ (7l |^n| 4 + 7|^n T) + f IV'n T |Cf ” ^ (|^n T ~ |C| 2 )} +( d - i ) 

for the inter-SOC system (7), where c.c. stands for the complex-conjugate expression. 

2.2. Stationary solutions 

Stationary solutions for the above systems are sought for as 00 = A n exp (— ipt) and 
0“ = R n exp (—ipt), where p is the chemical potential. The set of stationary equations 
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pAn = — C(A n+1 + A n _ i) — in(A n+ i — A n _i) + (o / 1 1 1 2 + fd\B n \ 2 )A n — flB n , 

pB n = — C(B n . |_i + B n _ i) + in(B n+ i — B n _ i) + (/3|A n | 2 + r )\B n \~)B n — D^k?) 

and its counterpart following from Eq. (7) is 

f^An = — C(A n+ 1 + A n _ i) + in{B n+ \ — B n _ i) + ( 7 i|A n | 2 + (3\B n \ 2 )A n — TlA n , 

f-Bn = — C(B n+ 1 + B n - 1 ) + in(A n+ 1 — A n _i) + (/3|A n |' + ^\B n \ 2 )B n + Df(?Ri) 

Our aim is to find stable discrete solitons for both systems, with the intra- and inter- 
SOC. To this end, stationary equations (12) and (13) were solved by means of a 
numerical algorithm based on the modified Powell minimization method [46]. Stability 
of so found discrete solitons was checked, in the framework of the linear stability 
analysis, by numerically solving the corresponding eigenvalue problem for modes of 
small perturbations. Finally, Eqs. ( 6 ) and (7) were directly simulated by dint of the 
Runge-Kutta procedure of the sixth order, [46], to verify stability properties predicted 
by the linear analysis. This was done by performing direct simulations with generic 
small random perturbations added to the initial conditions. 

2.3. The dispersion relation 

The search for localized solutions - here, complexes formed by discrete solitons in each 
BEC component - starts from the derivation of the underlying dispersion relations. To 
this end, we substitute = i/j + exp ( i(kn — pt)) and iJj~ = -0 _ exp {i{kn — pt)), with 
wavenumber k, into the linearized version of Eqs. ( 6 ) and (7). Straightforward algebra 
gives rise to the following dispersion relation for both intra-SOC and inter-SOC systems: 

P 12 = —2 C cos k ± \A~2 2 + 4k 2 sin 2 k. (14) 

Further analysis of the dispersion relation is worded in terms of the intra-SOC model, 
for which discrete solitons are actually produced below. 

In the absence of the linear coupling between the two components, i.e., 0 = 0, Eq. 

(14) reduces to a straightforward generalization of the commonly known spectrum of the 
discrete linear Schrodinger equation. In the present case, the spectrum corresponding 
to O = 0 contains the Bloch band, |/xi j2 1 < 2 \/C 2 + k 2 , sandwiched between two semi¬ 
infinite gaps. As follows from elementary analysis of Eq. (14), and shown in Fig. 1(a), 
the same remains qualitatively correct if the SOC is relatively weak, viz., at |0| < 2 C. 
Further increase of the strength of the linear coupling between the components opens 
a mini-gap in the middle of the Bloch band, which exists at |f2| > 2 C, occupying the 
region of 

M < M - 2C, (15) 

as shown in Fig. 1(b). Usually, mini-gaps open up in the spectrum of ID [47] and 2D 
[48] superlattices, i.e., discrete lattices which are subjected to an additional spatially 
periodic modulation. A mini-gap is also present in the spectrum of the discrete coupler 
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introduced in Ref. [49]. In the present case, it may be understood, roughly, as a result 
of the precession of the pseudo-spin around the magnetic axis. Lastly, Fig. 1(c) displays 
the limit case of C = 0, when the usual inter-site hoppings are negligible in comparison 
with the SOC-induced linear inter-site coupling, while the semi-infinite bandgaps and 
the mini-gap are still present in the spectrum. 

The tight-binding approximation underlying the derivation of the discrete equations 
takes into account only the first Bloch band, while the higher-order ones are neglected. 
This means that the band is sandwiched, as said above, between two semi-infinite 
bandgaps [see Fig. 1(a)]. In the presence of the attractive nonlinearity, the existence of 
fundamental ( unstaggered ) discrete solitons is expected in the lower semi-infinite gap. 
Solitons of this type are chiefly considered in this work. On the other hand, staggered 
solitons, with 7r-shifted fields at adjacent sites, may exist in the mini-gap [48, 49, 50]. 
Basic results for the latter type of the discrete solitons are briefly reported below too. 

3. Results and discussion 

In the following, results are chiefly presented for the intra-SOC model, based on Eq. (6). 
Qualitatively similar conclusions can be obtained for the inter-SOC model based on (7), 
which is explained by the fact that lattice solitons tend to be similar in different systems, 
due to the strong effect of the discretization. Recently, both 2D and ID discrete models 
with the inter-SOC terms (the Rashba coupling) were considered in [30]; however, unlike 
the present analysis, which focuses on strongly discrete states, that work was dealing 
with broad discrete solitons corresponding to the quasi-continuum regime. 

The discrete SOC-BEC systems support different types of bright solitons, depending 
on the systems’ parameters. They are classified here as miscible and immiscible 
complexes (MC and IMC, respectively, see Fig. 2), formed of localized patterns in both 
components [42, 43]. We consider as miscible the localized complexes in which both 
stationary components (A n and B n ) are overlapping, see Fig. 2(a). On the other hand, 
immiscible states are characterized by separated shapes of the two components. The 
numerical results demonstrate that stable IMCs are mostly formed by one component 
featuring two maxima and the other component with a single maximum between these 
two, see Fig. 2(b). It is relevant to mention that broad (quasi-continual) discrete 
solitons found for the inter-SOC system in [30] also feature miscible and immiscible 
structures, which are called, in that context, “smooth” and “striped” ones (the latter 
pattern demonstrates a large number of alternating local peaks of the density of the two 
components). 

Our analysis has identified four basic types of the MC solutions, which may be 
built as follows (cf. Ref. [51]): in-phase on-site solitons (IPOS), Re (A n ) x Re (B n ) > 
0, counter-phase on-site solitons (CPOS), Re (A n ) x Re (B n ) < 0, in-phase inter¬ 
site solitons (IPIS), Re (A n ) x Re (B n ) > 0, and counter-phase inter-site solitons, 
Re (A n ) x Re (B n ) < 0 (CPIS), see Fig. 3. 

Considering the IMCs, various structures of the on- and inter-site types have been 
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found, which reflects strong sensitivity of the system to parameter values. We report 
here only those IMCs which may be stable, and mention scenarios of the transformation 
of unstable IMCs into MCs. In fact, only on-site IMCs are found to be stable in certain 
ranges of the parameter space. 

Miscible on-site complexes of both the in-phase and counter-phase types are found 
to be stable over large parts of their existence region in parameter space ( 7 , /3, k, fl, C). 
As mentioned above, they belong to the lower semi-infinite gap of the systems’ spectra, 
shown in Fig. 1. IPOS complexes lose the stability only in a limited parameter area; 
for example, in the case of C — 1 they are unstable at j3 > —1 and —2 < O < —0.5. 
In direct simulations, the instability in that area leads to spontaneous formation of 
breathing localized structures with a smaller norm, in comparison with the original 
unstable discrete solitons. Without the inter-species attraction, /? = 0, IPOS solutions 
are unstable, except when the Rabi coupling is absent too, 0 = 0, making the system 
decoupled. On the other hand, CPOS are unstable at (3 < — 1 and 0^0. Without the 
attraction between the components, /3 = 0, and for small values of the Rabi coupling 
parameter (—1 < O < 0), CPOS are stable. For higher values of |0|, CPOS solutions 
break into breathing structures. Both components of the complex exhibit the same 
behavior and evolve into miscible breathers, exchanging the power (norm) in the course 
of the evolution. All other miscible solutions were found to be unstable. 

All inter-site-centered solitons are found to be unstable, as it usually happens in 
other discrete systems. A trend of the evolution of these unstable modes revealed by the 
simulations is to end up as on-site breathers with significant emission of matter waves. 
In particular, the CPIS complexes evolve so that both components exhibit the same 
behavior and evolve into miscible on-site breathing structures exchanging the power in 
the course of the evolution. In this case, a part of the emitted matter wave propagates 
in the form of a new two-component moving complex, see Fig. 4. 

Generally, we conclude that all unstable soliton solutions tune their energy, by 
radiating away a part of it, to the energy corresponding to the miscible 011 -site solitons 
existing at the same parameters . The modes evolve as stable on-site breathers. If the 
radiated energy, transferred to the background, is large enough, as is the case for the 
CPIS mentioned in the previous paragraph, new moving localized breathing structures 
can be formed (Fig. 4). 

The IMC solutions tend to get destabilized with the increase of the strength of 
the linear coupling between the components. Usually, they too evolve into miscible 
breathing modes after radiating away a part of the energy. Stable IMCs feature on-site 
two-component configurations, see Fig. 2(b). In the area of the system’s parameters 
where the IMCs are stable in the inter-SOC systems, a stable branch of MCs of the on¬ 
site type can be found too, which implies bistability between the IMC and MC modes. 
Then, selection of one of the coexisting solutions is determined by initial conditions. By 
checking values of the Hamiltonian, we have found that, in the bistability region (the 
whole gray area in Fig. 6 (b)), the MCs feature smaller energy, hence they realize the 
ground state of the system. The energy difference between IMC and MC in this case 
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is smaller for smaller values of C. Similar bistability is found in the inter-SOC system 
for the parameter set shown in Fig. 6(a) in a narrow area with very small « 0 and 

0 < C < 0.2. 

Unlike the intra-SOC system, in its inter-SOC counterpart the SOC can induce 
destabilization of the IMCs (leading to a spontaneous transition from the IMC to an 
MC) if one sets (3 = = C = 0 in Eq. (7), so that the two components interact 

solely via the linear SOC terms ~ k. This happens in a narrow area of the remaining 
parameter space. 

On the other hand, the SOC affects the threshold of the immiscibility-miscibility 
transition in both the intra-SOC and inter-SOC systems, in cases of the nonlinear 
coupling only (0 = 0), linear coupling only ((3 = 0), and when both couplings are 
present (O ^ 0 ,f3 0). Figure 6 displays the position of the miscibility threshold, 

i.e., the curve separating the region of the coexistence of stable IMC and MC, and 
the region where the IMC is unstable while the MC is stable (as concerns the on-site 
configurations), in the (|0|,C) parameter plane. In the inter-SOC system [Fig. 6(a)], 
parameter (3 = — 1 and k — 1 are fixed. In the intra-SOC system, we set f3 = 0 and 
k — 1, because for (3 < —0.5 stable IMC have not been found [Fig. 6(b)]. 

Figure 6 demonstrates that the IMC complexes have the largest stability region 
in the case of C = 0, when the usual hopping is negligible, while the SOC terms are 
present [recall that, as shown in Fig. 1(c), the respective linear spectrum remains 
meaningful, always featuring the mini-gap]. This case corresponds to the setting when 
the OL potential is very deep, and, simultaneously, the SOC, induced by the laser 
beams illuminating the condensate, is strong. Such a setting is feasible, because, in the 
underlying tight-binding approximation, the hopping coefficient, C in Eqs. (6) and (7), 
is proportional to the square of the small overlap integral of wave functions (Wannier 
modes [44]) localized at adjacent sites of the OL, while the scaled SOC coefficient, k, 
is proportional to the first power of the same integral. In fact, this limit case emulates 
the spin-orbit-coupled binary condensate made of infinitely heavy atoms (as mentioned 
above, the case of an infinite particle’s mass was recently considered in [15], in the 
context of a continual model elaborated for the measurement of spin). 

Furthermore, in the same limit of C —> 0, the staggering transformation in the form 
of if n = i n fi n [49], cf. Eq. (8), can be applied to the intra-SOC system, transforming 
it into a discrete system without the SOC terms, but, effectively, with the hopping 
coefficients of opposite signs, k and — k, in the two components: 

1 + V’n-l) + (Tl IV’n" I” + P$n\ 2 )$n ~ ^n , 

i ~W = ~ K ^ +l + ^-l) + Wn! 2 + 7l^n| 2 )^n ~ • (16) 

In this way, the model can emulate a two-component BEC with opposite signs of the 
effective atomic masses , which is not possible in the implementation of the BEC. 

Continuing the consideration of IMC modes in the intra-SOC system, we outline 
typical instability-development scenarios for them (similar scenarios are produced by the 
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inter-SOC system). Unstable on-site IMCs may actually be very robust overall, evolving 
into localized breathing miscible structures, see Fig. 7. In that case, the two components 
of the on-site IMC evolve into miscible breathing structures, which exchange the norm 
in the course of the evolution. 

Figure 8 illustrates the dynamics of an unstable inter-site IMC. In the course of 
the evolution, a large share of the norm is radiated away, and the formation of irregular 
localized patterns is observed in a neighborhood of the initial position of the components. 

All the localized complexes considered above are generated in the semi-infinite gap, 
in terms of the spectra of the linearized systems. It remains to briefly describe what 
happens when a mini-gap opens at |Q| > 2C, as shown in Figs. 1(b) and (c). In 
that gap, we have found miscible CPOS and CPIS complexes, which resemble similar 
complexes found in the semi-infinite gap, see Fig. 3. Only the CPOS complexes may 
be stable in certain parameter areas inside the mini-gap. All other structures are highly 
unstable in this case. 

To summarize this section, the stable on-site localized MCs (miscible complexes) 
can be created in the large area of the parameter space in the two-component SOC- 
BEC system. On the other hand, stable IMCs (immiscible complexes) were found only 
in a narrow parametric area. The latter one is often a bistability region, where the 
system supports both the IMC and MC stable complexes (which is usually not possible 
in the free space [30]). Unstable IMCs, as well as unstable inter-site MCs, evolve into 
breathing miscible states, emitting matter waves away. These conclusions pertain to 
both the intra-SOC and inter-SOC systems considered. Thus, we conclude that, in 
general, all inter-species couplings, linear (SOC) and attractive nonlinear ones, favor 
miscible states. 

4. Conclusions 

We have introduced a discrete one-dimensional model of SOC (spin-orbit-coupled) 
binary BEC trapped in a deep OL potential. Two systems were considered, of the intra- 
SOC and inter-SOC types, in which the SOC terms, represented by the discrete version 
of the first spatial derivatives, act, respectively, inside of each component corresponding 
to the dressed atomic state, or couple the different components. Each system includes 
attractive intra- and inter-species cubic interactions (repulsive interactions can be 
transformed into the same form by means of staggering). It has been demonstrated 
that both systems support two distinct types of discrete soliton complexes, miscible 
and immiscible ones, depending on values of the systems’ parameters. The transition 
between the two types of the complexes can be controlled by the SOC strength. As 
usual, only on-site-centered modes may be stable. The SOC terms may open up a 
mini-gap in the systems’ spectrum, where stable miscible on-site soliton complexes exist 
too. 

New possibilities offered by these settings are scenarios for the emulation of SOC 
binary condensates built of infinitely heavy atoms, as well as of the binary BEC with 
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effective atomic masses which have opposite signs. Both settings, that correspond to the 
limit when the usual hoppings are negligible, while the SOC keeps acting, are obviously 
unavailable in the direct realization of BEC. 

To continue the present analysis, it may be interesting to consider bound states of 
the discrete solitons, as well as their mobility. A challenging problem is to extend the 
analysis to the 2D setting, where it has been recently demonstrated that the SOC, in 
the combination with the nonlinear attraction, supports stable bright solitons in the in 
the presence of the OL potential [27, 30], as well as in free space [28, 30] (without the 
SOC, all 2D solitons sustained by the cubic attractive interaction in the free space are 
unstable against the critical collapse). 
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Figure 1. (Color online) Dispersion curves, as given by Eq. (14), for (a) ft = —0.5, 
C = 1; (b) ft = —6, C = 1; and (c) ft = —0.5, C = 0. In all the plots, k = 1. At 
|ft| > 2C, the mini-gap Eq. (15) opens up around p = 0. This figure corresponds 
to both systems considered here, of the intra-SOC and inter-SOC types, i.e., given by 
Eqs. (6) and (7), respectively. 
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Figure 2. (Color online) Panels (a) and (b) displaye typical shapes of the discrete 
miscible and immiscible complexes (MCs and IMCs), respectively. Black (squares) and 
red (circles) lines denote, severally, components A n and B n . The respective parameters 
in Eq. (12) are re = 1, 71 = —1, 7 = —0.5, /? = —1, f1 = —2, p = —19.4, C = 1 . 
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Figure 3. (Color online) Examples of discrete two-component solitons of the miscible 
type: (a) an IPOS complex, (b) a CPOS complex, (c) an IPIS complex, and (d) a 
CPIS complex. Black (squares), red (circles), and blue (triangles) lines depict real and 
imaginary parts and the absolute value of the fields, respectively. The parameters are 
k = 1, 7 = 7 i = — 1 , C = 1, fi = —9.4. 
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Figure 4. The evolution of two components of an inter-site CPIS mode with p = —9.4, 
n = 1, 7 = —0.5, 71 = —1, /? = —2,12 = —0.5, and C = 1. Both components break 
into a trapped on-site mode and a moving one. Initial small random perturbations are 
added to the solutions to catalyze the onset of the instability. 




Figure 5. (Color online) An example of a stable IMC (immiscible complex): (a) real 
parts of its two components; (b) the corresponding absolute values. The parameters 
are: II = —4, 71 = —1, 7 = —0.5, /? = —5, p = —16, k = 1 and (7=1. Black 
lines with squares and red lines with circles are associated with components A and B 1 
respectively. 
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Figure 6. (Color online) The transition from stable (gray areas) to unstable (white 
areas) IMC in the (101,(7) plane: (a) in the inter-SOC system [Eq. (7) with /3 = — 1], 
and (b) in the intra-SOC one [Eq. ( 6 ), with (3 = 0]. Other parameters are 7 = 71 = —1, 
and k = 1. In the area with small |0| in the grey area in plot (a) stable IMC and MC 
can exist. While in the plot (b) in the entire gray area both stable IMCs and MCs 
can be created (bistability regions). Those IMCs which are unstable (white areas in 
figure) spontaneously evolve into stable MC breathers, after radiating away a part of 
the initial norm. 




n n 

Figure 7. The evolution of an unstable immiscible complex (IMC) in the intra-SOC 
system with // = —19.4, n = 1, 7 = —0.5, 71 = —1, /3 = —1,0 = —2, and (7=1. 
Initial small random perturbations were added to the immiscible localized complex. 
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Figure 8. (Color online) The evolution of an unstable inter-site IMC in the intra-SOC 
system. The parameters are the same as in Fig. 7. 






